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We study the distance dependence of the rate of electronic excitation energy trans- 
fer from a dye molecule to a metal nanoparticle. Using the spherical jellium model, 
we evaluate the rates corresponding to the excitation of I = 1, 2, and 3 modes of the 
nanoparticle. Our calculation takes into account both the electron-hole pair and the 
plasmon excitations of the nanoparticle. The rate follows conventional i? -6 depen- 
dence at large distances while small deviations from this behavior are observed at 
shorter distances. Within the framework of the jellium model, it is not possible to 
attribute the experimentally observed d~ 4 dependence of the rate to energy transfer 
to plasmons or e-h pair excitations. 

I. INTRODUCTION 



Fluorescence Resonance Energy Transfer (FRET) is an interesting photo-physical pro- 
cess l) that involves transfer of excitation energy from donor to acceptor in a non-radiative 
fashion. It has been extensively used in biology as a spectroscopic ruler to study the confor- 
mational dynamics of biopolymers in the 10-100 A range. The rate of non-radiative damping 
of dye molecules near quencher molecules is found to vary as R~ 6 , where R is the distance 
between the donor and the acceptor. In FRET, energy is transferred from a donor molecule 
to an acceptor molecule via the dipole-dipole interaction. The matrix element for the in- 
teraction is given by Vda oc ^ D ^lf a where Jad and /2a are the transition dipole moments 
of donor and acceptor. Therefore, the rate of energy transfer, which is proportional to the 
square of the interaction matrix element varies as R~ 6 . The Forster expression for the rate 
is given by k nr = k r (j^j where k r is the radiative rate and R is the Forster radius, which 
can be expressed as an overlap integral between the donor emission and acceptor absorp- 
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tion spectra as R oc (^k 2 / ^1 



d (v) a A W, 



The orientation factor k 2 takes into account 



the effect of relative orientation of the donor and acceptor transition dipole moments and 
is given by k 2 = (sin 6 A sin 6b cos {<pa — <Pd) — 2 cos 6a cos 6d) 2 . The value of k 2 is usually 
taken to be | assuming rapid orientational averaging of the donor within the lifetime of its 
excited state Q]. 

However, there are reports of deviations from the conventional FRET behaviour in the 
literature. Two possible reasons for this are the breakdown of point drpole approximation, 
especially for systems with extended transition densities like polymers [3( and the incomplete 
orientational averaging within the lifetime of the excited state 4|. Currently, there is a 
lot of interest in using nanoparticles [s| as quenching agents due to their tunable optical 
properties 6j. There have been very interesting observations on energy transfer between the 
dye fluorescein and a gold nanoparticle of diameter 1.4 nm. Using double stranded DNA as 
the spacer, the dye molecule was kept at different fixed distances from the particle and the 
rate of energy transfer was experimentally determined [l, 8]. They found a d~ 4 dependence 
of the rate on the distance d between the dye and the surface of the nanoparticle and have 
referred to this as nanoparticle surface energy transfer (NSET). Such a dependence is of 
great interest as it would more than double the range of distances that can be measured. 
Persson and Lang (PL) |9| had long ago studied the dissipation of vibrational energy of an 
oscillating dipole held at a distance d above the surface of a semi-infinite metal. In this case, 
the energy is transferred to electron-hole (e-h) pair excitations in the metal, which form 
a continuum, having all possible energies from zero to infinity. For an excitation energy 
u, their density of^states is proportional to uj. PL found the distance dependence to be 
H]). Following this, it is suggested 3, s| that the observed distance 



d 4 (see also 
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dependence in NSET is due to the excitation of e-h pairs in the nanoparticle. In support of 



this, the authors [7|, [8J point out that plasmonic absorptions have not been observed for a 
gold nanoparticle of diameter 1.4 nm. However, unlike a semi-infinite metal, a nanoparticle 
of such a small size, does not have a continuum of e-h excitations. The excitations are 
discrete (see the excitation spectrum given in Fig. [2] of this paper). In view of this, it would 
be very interesting to study the distance dependence of NSET theoretically. 

There have been attempts to explain the d~ A dependence for the case of nanoparticle 
[12] . [lI \l4 \. In particular, focus of these papers [ljjl is on energy transfer to the plasmons of 
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the nanoparticle and they find a predominantly R~ 6 dependence, though deviations occur 
at shorter distances. These calculations do not account for e-h pair excitations. To the best 
of our knowledge, there are no reports of theoretical calculations which take plasmons as 
well as e-h excitations of the nanoparticle into account. Therefore, we have calculated the 
rate of transfer of electronic excitation from an excited fluorescein molecule to a 1.4 nm 
nanoparticle, taking both the e-h pair excitations as well as the plasmon excitations of the 



nanoparticle into account. We model the nanoparticle within the jellium model [15], and 
use time dependent local density approximation (TDLDA) to calculate the excitations of 
the nanoparticle. Strictly speaking, one should use RPA with the Hartree Fock approach 
as the starting point, because in any LDA based approach, the single particle energies and 
the response calculated using them have no fundamental meaning. However, there have 



been a number of calcu 



ations based on TDLDA for the response of a metallic nanoparticle 



and other systems 0, 0, Q, Q- 

All these calculations have led to results which agree 
well with experiments and therefore we use this approach. Within this framework, we do 
not get the experimentally observed d~ A dependence. Therefore, in our opinion, this very 
interesting experimental observation is, as it stands, unexplained. Perhaps, it may be due to 
other factors, like transfer through the DNA, or the asphericity of the nanoparticle, which 
are being currently investigated. It is also to be noted that the distances between the donor 
and the acceptor in the NSET experiment are much greater than the dimensions of the 
donor and acceptor, making it unlikely that the breakdown of point dipole approximation 
is responsible for the d~ 4 dependence. 

II. MODEL FOR THE NANOPARTICLE 

We use the spherical jellium model for the nanoparticle which provides a model system 
for investigating the response of the conduction electrons in small metal particles 2(3, 21]. In 



this model, the positive ions of the metal cluster are replaced by a uniform sphere of positive 
charge and the density functional formalism within the linear response approximation is used 
to calculate the response of the cluster to the time dependent external potential. Two kinds 
of excitations are possible for a metal cluster namely, single particle excitations and plasmon 
excitations. In a single particle excitation, an electron is excited from an occupied level to 
an unoccupied level and these are the electron-hole pair excitations. Plasmon excitations 
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are collective oscillations involving many electrons, wherein the electronic charge density 
oscillates as a whole against the positive background. 

We denote the Hamiltonian of the nanoparticle as Hq. The molecule is treated within 
a single particle model. An electron which was initially in the orbital <fi g is excited to an 
orbital <fi e . De-excitation of the molecule, in which the electron goes back to <p g may be 
thought of as a time dependent potential acting on the nanoparticle, given by 

w)-^/«gf*gV (i) 

J \r — r | 

In the above, u is the frequency of the transition corresponding to de-excitation. Thus the 
total Hamiltonian can be written as: 

H = H + <&ext{r, t). (2) 
One can carry out an expansion of the electrostatic potential at f due to a charge distribution 
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t—7] = E ^-mYUmAti), (3) 

I I l,m > 

where r> (r<) is the larger (smaller) of r and r and Y^ m (Vt) are the spherical harmonics. 
From Fig. [H it is obvious that the integration over r is to be performed over the molecule. 
As the electron density of the molecule is fully outside the nanoparticle, r < r and hence 

$ ex t(r,t) = e-^Y.M^Y*^) (4) 

l,m 

where 



4?r /■ <j}* e {r') (j)g{r) 
21 + 1 J r" 

Thus the perturbation acting on the nanoparticle is a combination of multipole fields of 
the form r l Yf m {VL) for various values of I and m. The case of such perturbations acting on 



nanospheres of jellium is well studied in the literature 
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III. THE RATE OF ENERGY TRANSFER 



The rate of transfer using the Fermi golden rule is: 

k = Y E I ( E \®™t{r, t)\G) \ 2 5(E E -E G - hu) (6) 

E 
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where \G) is the initial (ground) and \E) is the final (excited) state of the nanoparticle. For 



a spherically symmetric closed-shell system 



24|, within the jellium model, the perturbation 



r l Y* m (p) can only lead to the excitation of the lm th mode in the nanoparticle. Thus, the 
excitations in the nanoparticle have the same symmetry as the perturbation and the rate 
simplifies to 

k = ¥ E M U (E^Y^mG) \ 2 5(E E ^ -E G - hu). (7) 

11 Elm 

\Ei >m ) denotes excited states having quantum numbers l,m. 

If the nanoparticle is placed in an oscillatory external field & ext (r)e~ luJt = r l e~ lwt , an 
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induced electronic charge density is set up, which within the linear response theory 
is given by 

5 P i{r) = J dr'Ui(r,r\oj)^ l ext (r') (8) 

where II; (r, r\uj) is the density-density correlation function or the polarization propagator. 
Within the random phase approximation (RPA), 11/ (r, r',u) obeys the following integral 



equation: 
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Ui(r,r',uj) = ltf(r,r',uj) + J dr 1 J dr 2 n°(r, r u u)Ki{r x , r 2 )Ui(r2, r, u) (9) 

where n°(r, r , to) is the independent particle propagator and /^(ri,r2) is the effective two 
particle interaction [2l|. The independent particle approximation to the response function 
Il^r, r , lu), contains the e-h pair excitations only, while the response under RPA, given by 
II;(r, r ,uS) includes both the single particle and plasmonic response. The free and RPA 
response to the external perturbation of the system is calculated by integrating the external 



potential over the transition density 
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21] . Thus we have 



n°«,,^)= J dr$ l ext (r)SpKr) (10) 

and 

nf PA (<&U w) = / dr* l ext (r)Sp? P A (r) (11) 



The polarization propagator is related to the strength function, Si(u) by 

S^uj) = -Imll^i^u). (12) 

7T 



The rate may be expressed in terms of the strength function Si(u) as 

27, 



^tE m ;» ( i3 ) 

E,lm 
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The evaluation of polarization propagator for the independent particle as well as the RPA 
response allows one to take both the e-h pair excitations as well as the plasmonic excitations 
into account in the calculation. We make use of the time dependent local density approxi- 
mation (TDLDA) version of the above equations, as this has been found to lead to excellent 
results. Details of the approach may be found in 17|, |21|. 



IV. CALCULATIONS 



We have optimized the geometry of fluorescein, using the Gaussian03 program, within 
the DFT approximation (B3LYP-6-31G*). The HOMO (<p g ) and LUMO (0 e ) were taken 
from the calculation for the optimised geometry. These orbitals were then used to evaluate 
the matrix elements M/ m numerically for all m with I up to 3. A grid within a box of 
size 24A x 24A x 24A, within which the molecule was located was used for this purpose. 
Note that / = 1 corresponds to the oscillation of the electrons in the nanoparticle that 
has the shape of a p-orbital, I = 2 that of a d-orbital and I = 3 that of an f-orbital. We 
assumed an Au cluster of 90 atoms which corresponds to a 1.4 nm gold nanoparticle and 
performed the jellium model calculation to evaluate Itf(& ext , to) and ITf p {& ext , oj). The rate 
was calculated using Eq. (TlBl . with static but random averaging over all the orientations of 
the nanoparticle. Thus, we have evaluated the rate of energy transfer from the fluorescein 
molecule to the 1.4 nm gold nanoparticle as a function of distance. 



V. RESULTS AND DISCUSSION 



The rate of energy transfer depends on the values of the two terms M^ m and Si(u>). 
The first one represents the perturbing potential acting on the nanoparticle due to the dye 
molecule and the second one is the response of the nanoparticle to the perturbation. We 



have evaluated SAu) for I — 1, 2 and 3 modes of excitations o 



TDLDA response using the JELLYRPA program of Bertsch 21| . For the calculation, we 



the nanoparticle, for the 



have taken r s = 3 atomic units (au). The broadening parameter T, which determines the 
width of the single particle peaks was taken to be 0.01 eV. Plot of Si(u) calculated using the 
single particle and RPA approaches is shown in Fig. [2]for I = 1, 2 and 3 modes. The TDLDA 
plot clearly shows peaks corresponding to e-h pair excitations, collective surface plasmon 
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SAlo) (SP) 


5,(a;) (TDLDA) 




(a.u. 2 7eF) 


(a.u. 2 7ey) 


1 


32.69 


14.32 


2 


259500. 


322.4 


3 


5.714 x 10 6 


3.988 x 10 6 



TABLE I: Numerical values of Si(u) for / =1, 2 and 3 at w = 2AeV, the emission energy for 
fluorescein. SP stands for single particle. 

and bulk plasmon excitations. The numerous narrow spikes at lower energies correspond to 
e-h pair excitations. The peak around uo/ujb — 0.5 is the remnant of the surface plasmon and 
the one around oj/ujb — 1-2 corresponds to the bulk plasmon. Note that the frequency of a 



surface plasmon mode |25j is related to that of the bulk plasmon mode as uo s = <jJb 
The surface plasmon mode is slightly red shifted while the bulk plasmon mode is sli 
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blue shifted. The shifts are consistent with previous jellium model calculations 
The quantities of interest to us, in the calculation of the rate of energy transfer to the 
nanoparticle are Si(u) for I = 1, 2 and 3, with u having the value of emission energy of the 
dye particle, which are given in Table 1 for uj = 2AeV . 



After evaluating M^ m numerically for all m with I upto 3, we evaluated the rates of transfer 
corresponding to the I = 1, 2 and 3 modes of excitation of the nanoparticle, both for the 
single particle as well as the TDLDA response. The results for TDLDA are shown in Fig. 
[3j Note that the expansion in Eq. ([3]) and hence the calculations are not valid if the dye 
penetrates into the electron density of the nanoparticle. Moreover, in such a case, overlap 
effects dominate and one has to think of the Dexter mechanism and evaluate the exchange 
integral. The rates for I = 1, 2 and 3 modes vary with distance as i?~ 6 , R~ 8 and R~ 10 as 
may be seen from the slopes in Fig. [31 Then, we evaluated the total rate of transfer, again 
for single particle as well as TDLDA response (see Fig. H]). It is found from the figure that 
the single particle rate is actually larger than the TDLDA result. This is because there 
is intensity borrowing by the collective modes (bulk and surface plasmons) from the single 
particle excitations, resulting in a lower value for Si(u) in the TDLDA calculations, as is 
clear from table fl] Asymptotically, at large distances, it is only the / = 1 mode that is 
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important and this leads to the conventional ~ R~ G limit for the total rate. But, at shorter 
distances, I = 2 and 3 modes gain importance leading to deviations from R~ 6 behavior at 
shorter distances. To make this clearer, we fitted the long distance rate with and used 
the result to calculate the short distance rate. The resultant rate is shown in Fig. HJ along 
with the actual rate. The actual rates are found to be slightly higher than expected from 
~ R~ 6 dependence. This is due to the fact that I = 2 and 3 modes become important at 
lower distances. Note that the deviations are not such as to give an R~ A dependence in the 
distance range that we study. It is also of interest to note that the R dependence is governed 
by |M; jm | 2 and not by Si(u). Use of a different set of values of r s , or uo would not change 
the R dependence. Therefore, we conclude that excitation of plasmons or e-h pairs cannot 
lead to the observed experimental data. In ref. 14j], for a distance range from infinity to 
up to 4 times the radius of the nanoparticle (for our case, 28 A), the behavior is found to 
be R~ G . In this range, we too get an R 6 dependence, illustrating that inclusion of e-h pair 
excitations does not affect the result at all. Also, it is to be noted that in our calculations, 
our grid used for numerical integration will penetrate into the nanoparticle for distances less 
than about 28 A so that our approach will break down for closer distances. 



VI. CONCLUSIONS 



We have adopted a spherical jellium model to evaluate both the independent particle and 
collective response of a metallic nanoparticle to an external time dependent perturbation 
and used this to evaluate the rates of non-radiative energy transfer from the excited state of 
a fluorescent dye to the particle. The rates due to the excitation of I = 1, 2 and 3 modes of 
the nanoparticle were found to vary with distances as R~ 6 , R~ 8 and R~ 10 respectively. The 
major contribution to the rate is from the leading R~ 6 term. The contributions from the 
other terms are relatively small, except at distances close to 28A, and these make the rate 
only larger and they are not in a direction as to make the result behave like d~ A . Thus the 
experimentally observed d~ A dependence cannot be explained by considering the excitation 
of plasmons or e-h pairs of the nanoparticle. 
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FIG. 1: A schematic of the system consisting of the gold nanoparticle and the dye molecule. 




FIG. 2: The imaginary part of the RPA polarization propagator for the excitation of I =1, 2 and 
3 modes of the nanoparticle. The frequency is given in units of the bulk plasmon frequency. 
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FIG. 3: Distance dependence of the rate of transfer corresponding to the excitation of / =1, 2 and 
3 modes of the nanoparticle for the TDLDA response. The slopes of the log-log plots are —6.005, 
-8.019 and -10.039 for I =1, 2 and 3 respectively. 
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FIG. 4: Distance dependence of the calculated total rate of transfer for the single particle and 
RPA responses. In both the cases, the rate calculated using asymptotic expression c/R 6 (referred 
to as asymptotic rates) are shown as dotted lines. The actual rate is greater than this due to the 
contributions from I = 2 and I = 3 modes. 



